Spin bags in the doped t—J model 
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We present a nonperturbative method for deriving a quasiparticle description of 
the low-energy excitations in the t—J model for strongly correlated electrons. Using 
the exact diagonalization technique we evaluated exactly the spectral functions of 
composite operators, which describe an electron or hole dressed by antiferromagnetic 
spin fluctuations as expected in the string or spin bag picture. For hole doping up 
to 1/8, use of the composite operators leads to a drastic simplification of the single 
particle spectral function: at half-filling it takes free-particle form, for the doped case 
it resembles a system of weakly interacting Fermions corresponding to the doped 
holes. We conclude that for all doping levels under study, the elementary electronic 
excitations next to the Fermi level are adequately described by the antiferromagnetic 
spin fluctuation picture. 
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Despite great efforts the unusual properties of high-temperature superconductors remain 
a largely unresolved problem. There is general agreement that the electrons in these mate- 
rials are strongly correlated, so that marked deviations from the single particle picture are 
to be expected. An important step in setting up a successfull theory of both, their normal 
and superconducting state, therefore would be to find a description in terms of elementary 
excitations i.e. weakly or moderately interacting quasiparticles. We present an exact diago- 
nalization study of the t—J model which is specifically aimed at finding such a description. 
We construct composite operators which reduce or increase the electron number by one and 
simultaneously rearrange the spins in the neighborhood of the newly created hole/electron 
so as to simulate the 'cloud of spin defects' surrounding the hole. For doping levels up 
to 1/8, the spectral function of the composite operators then takes the form expected for 
weakly interacting Fermions, with a well defined 'band' right at the Fermi level. The t—J 
model reads: 

H = -t Y: (4.^,,. + H.c.) + JY.\S^-S,-'^]. 

<i,j>,c7 <i,j> 

The Si are the electronic spin operators and the sum over < i, j > stands for a summation 
over all pairs of nearest neighbors on a two dimensional square lattice. The operators cl^. 
are expressed in terms of ordinary fermion operators as c\ ^{1 — 71^ _^). We present results for 
a 4x4 cluster of this model with t/J=A, similar results have been obtained also for different 
values of t/ J. 

A single hole in the half-filled band can be well described by the string [||-|^ or spin bag 0] 
picture, where the hole is dressed by antiferromagnetic spin fluctuations. With this in mind, 
we make the following ansatz for a 'spin bag operator' {N{j) denotes the nearest neighbors 
ofj): 

= E ^i^k,i, 
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keN(j) l£N{k) 

A-,3= E E E S^S^SfCn^^i. (1) 
keN{j) l£N{k) meN(l) 

When acting on the Neel state, the operator Aj x creates all strings of length A which begin 
at site j, Ck,u reproduces a simple trial wave function for a single hole p. Aj^x also can be 
thought of as having been generated by A-fold commutation of Cj^-^ with the kinetic energy, 
a procedure suggested by Dagotto and Schrieffer |Q. The parameters ax{k) are determined 
from the requirement that the state Cfc_o- |^o°'^^) (where I^I'o'"*'') denotes the half- filled 
ground state) has norm 1 and maximum overlap with the lowest totally symmetric single- 
hole eigenstate with momentum fc, |\E'o^''^(fc)). If we denote the Fourier transform of Aj^j, by 
Ah y and introduce 

n. = (vl/(^'^)(fc)|A,,.|vl/r)), 



5M,.=(^ri4,,^Mi^r), (2) 

the a's can be obtained by solution of the generalized eigenvalue problem Ax = XBx. We 
use Amaa;=3, SO that the Cfc,o- effectively contain 3 free parameters. Having fixed the a's, we 
compute the spectral function 

A.{k, -co) = T.\{^nc.,.\<'^)mu; - {El^^^ - E^^)), 

V 

via the standard Lanczos method. Here \^^^^^) denotes the v^^ eigenstate with n holes and 
^(nft.) corresponding energy (in particular z/ = implies the ground state). This spectral 
function is shown in Fig. |l| and compared to that of the 'string-0 operator' V2ck,a, which 
up to a factor 2 equals the usual photoemission spectrum. Despite the small number of 
adjustable parameters, use of the spin bag operators brings the spectral function to almost 
perfect free-particle form: the incoherent continua present in the bare electron spectra are 
removed, the spectral weight being concentrated essentially in one sharp peak with a well- 
defined next-nearest neighbor dispersion. This can be understood if one attributes the 



incoherent continua to the retracable-path-type motion |^ of the hole. Since the Aj^^ describe 
precisely the forward and backward hopping of the hole along a track of Neel ordered spins, 
this type of hole-motion is already 'incorporated' into the definition of the Ck^a, so that their 
spectra emphasize the coherent hole motion. Some care is necessary: for momentum (vr, vr) 
the lowest totally symmetric single-hole eigenstate has spin |, so that it can not be observed 
in the spectrum of the bare electron operator. On the other hand, since the cg^ are not 
vector operators under spin rotations, it can be observed in their spectrum. Since this state 
fits very well into the next-nearest neighbor hopping dispersion relation expected for a single 
hole, we believe that the fact that it has spin 3/2 in the 4x4 cluster is a finite size effect. We 
do not have any proof for that, so results for (vr, vr) should be considered with care; however, 
none of the conclusions to be presented below depends crucially on the form of the spectra 
for this momentum. 

So far, we have merely demonstrated the quality of the string picture at half-filling which 
may not be very surprising; a much more interesting question is, whether this description of 
the states next to the Fermi level remains valid upon doping. Exact diagonalization offers 
a very direct and natural way to check this ssue, namely to evaluate the spectra of the ^ 
for a doped rather than the half-filled ground state. For the a\{k)^s we thereby retain the 
values optimized at half-filling (it would be easy to recalculate the a's for the doped ground 
state but an important question is whether there is some continuity in the development of 
the electronic states at the Fermi level). We begin with the ground state with 8 up-spin 
electrons, 7 down-spin electrons and momentum (7r/2, 7r/2)) i.e. the single hole ground state. 
In Fig. (Fig. H) the spectra of cg^ (cg j) are again compared to those of the respective 
string-0 operators, V^cj:^ (\/2cg^). Quite obviously, the c^^ continue to 'work': there is 
the same elimination of the incoherent continua and enhancement of the peaks at the Fermi 
level as for half-filling. A novel feature is the broadening of these peaks, and it seems natural 
to ascribe it to the scattering of the added hole from the one already present in the system. 
Most important of all, however, the spectrum for cj: ^ shows an unambiguous 'pocket' at 
(7r/2,7r/2), the momentum of the |-hole already present in the system: we clearly see the 
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Pauli principle working for the spin bags. This suggest weakly interacting spin- 1/2 Fermions 
which correspond to the doped holes as the 'effective theory' for the low-lying states of the 
cluster with 2 holes. 

We thus push things further and proceed to the two-hole ground state. Having in mind the 
results obtained so far, we should model it as an interacting state of two spin bags with total 
momentum zero: 

|§D = EA(fc)£M£-.al^r)- (3) 

k 

Using the Cfc,o- optimized at half- filling we consequently construct the states = 
Cfc,tc_fc4|^E'o°^^) for all 16 allowed momenta in the cluster, evaluate the matrices h^.y = 
{^{k)\H\^{k')) and nfcfc/ = ($(fc)| and solve the resulting eigenvalue problem to 

obtain A(fc). The estimate for the ground state energy obtained in this way is —8.21t, to 
be compared with the exact value of — 8.81t, and we have {^'^'^''1'$'^^^) = 0.67. Obviously 
the string gound state is not really an excellent approximation to the exact one, but it 
should be noted that it has been constructed from the half-filled groundstate, so that no 
relaxation of the 'spin background' (corresponding to the collapse of long range Neel order 
in the infinite system) is incorporated. On the other hand, the approximate ground state 
shares some basic features of the exact one, such as the correct (i2.2_y2-symmetry ^j. In Tab. 
|, the 'spin bag momentum distribution' n{k) = |($(fe)|^Q^''^)p/($(fc)|$(fc)), is listed; due 
to the nonorthogonality of the \^{k)) there exists no simple sum rule for this quantity, so 
that its interpretation as 'momentum distribution' is strictly speaking questionable. How- 
ever, it can give a rough idea of the distribution of the spin bags in momentum space, and 
obviuosly only (vr, 0) and (7r,7r/2) are appreciably occupied 0. It is interesting to contrast 
h{k) with the 'bare hole momentum distribution' nh{k) = {cj^^cl ) also given in Tab. | 
both for the approximate and exact two- hole ground state. Whereas the nh{k) for the exact 
and approximate ground state agree reasonably well, there is no similarity with h{k): nh{k) 
is roughly consistent with a free electron picture, n{k) would rather suggest that mainly 
(vr, 0) is occupied by quasiparticles. This can be understood by recalling that the incoherent 
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retraceable path- type motion of the bare holes (which naturally contributes to nh{k) 0) is 
already absorbed into the definition of the string operators, so that n{k) measures predom- 
inantly the coherent hole motion. 

Next, let us find an approximate annihilation operator for the dressed holes, i.e. an operator 
cl. ^ so that this operator and Ck^„ mutually undo the action of each other. Thereby the first 
guess, cl^a^ick^aV Hiay not be expected to be reasonable because the 'basis set' Aj^\ does not 
consist of operators which create electrons in orthogonal basis functions, i.e. [Aj,^, (A^ + 
We thus introduce a shortcut: for ^ we make the ansatz 

and construct the state 

If ^ and Cfc,o- indeed were adjoint fermionic creation and annihilation operators, the anti- 
commutator on the right hand side would be 1 and consequently we determine the /3's from 
the requirement that I^E'^^^-') has norm 1 and maximum overlap with |\['g^'^^). The actual 
value of the overlap then also provides a measure for the quality of our ansatz and with 
Amax = 3 we indeed find that for most fc-points (see Tab. 0), so that 

this way of obtaining an approximation to ^ appears quite reasonable. Using the ^ we 
now can also study the spin bag removal (electron addition) spectra in the two-hole ground 
state, defined as 

u 

Thereby some care is necessary: in contrast to the ordinary electron operators, there exists 
no simple sum rule for the integrated weight of the spin bag addition and removal spectra. 
To faciliate the comparison with the usual single particle spectral function in the following 
all spectra are therefore normalized to unity. 

Fig. ^ then compares the spin bag spectral function and that of the ordinary electron oper- 
ators. The spin bag creation (electron annihilation) spectra clearly show the simplification 
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already familiar from the previous calculations: the incoherent continua far from the Fermi 
energy are removed, the peaks near Ef are markedly enhanced. We conclude that even for 
this level of doping, the Cfc^o- optimized at half-filling provide a good description of the elec- 
tronic excitations closest to the Fermi energy, Ef. The degree of broadening of the spin-bag 
peaks is reminiscent of a Fermi liquid: there are sharp peaks close to Ef, diffuse peaks at 
lower energies. For the spin bag annihilation (electron creation) spectra use of the adjoint 
spin bag operators leads to an increase of spectral weight at (vr, 0) and a marked depletion 
at (7r,7r). The resulting division of spectral weight between spin bag addition and removal 
spectra moreover is remarkably consistent with the 'momentum distribution' n{k): spin 
bags can be annihilated predominantly at k={7T, 0) and k={TT, tt/2), i.e. the momenta which 
were most probable in the variational ground state I^q^'^''). All in all the doped cluster thus 
behaves very much like a system of weakly interacting 'effective fermions' corresponding to 
the doped holes. 

In summary we have shown that for all doping levels studied the use of 'spin bag operators' 
leads to a drastic simplification of the spectral function. Since these operators describe the 
modification of the 'spin background' in the immediate neighborhood of the hole, we also 
expect that they are predominantly determined by the local spin correlations, which most 
probably are described adequately by the exact diagonalization. It thus seems reasonable 
to expect that a similar simplification of the spectral function occurs also in the infinite 
system. Our nonperturbative approach thus suggests a rather conventional 'effective the- 
ory', in which the electronic excitations right at the Fermi level are modelled by spin 1/2 
Fermions corresponding to the doped holes (of course, like any other cluster calculation the 
present study cannot adress the dependence of the quasiparticle lifetime on the distance in 
energy from the Fermi energy). The remarkable degree of continuity upon doping moreover 
shows that these quasiparticles are holes dressed by antiferromagnetic spin fluctuations in 
essentially the same way as a single hole moving in an antiferromagnet, a reasonably well 
understood problem |^-^,[TD|JTT|| . The dressing of the holes with spin fluctuations leads to 
an interaction which favours a bound state with dx2_y2 symmetry. Our exact results thus 
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clearly corroborate the basic assumptions of the antiferromagnetic spin fluctuation theory of 
high-temperature superconductivity |T2| , p!3| . Another important point is, that the absorp- 
tion of the incoherent hole motion into the definition of the spin bag operators leads to an 
almost complete 'decoupling' of the quasiparticle spectral function and momentum distri- 
bution from that of the bare electrons. This suggests consequences for the interpretation 



of experiments: it seems plausible that measurements of transport properties |l4| , p!5| would 
rather probe the 'quasiparticle properties' and not resolve the internal structure of the spin 
bags. On the other hand, high-energy experiments like photoemission [l^Jl^l should resolve 
the internal structure and refiect properties of the bare electrons. Discrepancies between the 
transport properties and the Fermi surface measured in photoemission [|17| thus may be not 
surprising. 

It is a pleasure to acknowledge numerous instructive discussions with Professor S. Maekawa. 
Financial support by the Japanese Society for the Promotion of Science is most gratefully 
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FIGURES 

FIG. 1. Spectral function of the spin bag operators (full line) and ordinary electron annihilation 
operator (dotted line) in the half- filled ground state of the 4x4 cluster. 

FIG. 2. Spectral function A-{k,uj) of the spin bag operators (full line) and ordinary electron 
annihilation operator (dotted line) in the single hole ground state (momentum (7r/2, 7r/2)) of the 
4x4 cluster. The spin of the newly created hole is antiparallel to that of the hole already present. 

FIG. 3. Spectral function A_ (fe, a;) of the spin bag operators (full line) and ordinary electron 
annihilation operator (dotted line) in the single hole ground state (momentum (7r/2, 7r/2)) of the 
4x4 cluster. The spin of the newly created hole is parallel to that of the hole already present. 

FIG. 4. Spectral function A-{k, u) +A^{k, oj) of the spin bag operators (full line) and ordinary 
electron operators (dotted line) in the two hole ground state of the 4x4 cluster. The Fermi energy 
is marked by a thin line and the frequency region ooKEp (w>£^i?) corresponds to the annihilation 
(creation) of an electron. 
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TABLES 

TABLE I. 'Quasiparticle momentum distribution' h{k) in the approximate two-hole ground 

state and bare hole distribution function rihik) for the approximate (A) and exact (E) two-hole 
ground state. 



k 


(0,0) 


(i,o) 


(7r,0) 


/ 7r TT \ 

V 2 ' 2 / 




(7r,7r) 


n{k) 
nh{k),A 
nh{k),E 


0.0000 
0.0164 
0.0069 


0.0586 
0.0237 
0.0319 


0.7649 
0.2538 
0.1752 


0.0000 
0.0658 
0.0660 


0.2269 
0.2212 
0.2529 


0.0000 
0.2239 
0.2675 






TABLE II. The quantity 


'^^Q^^)^ for all momenta in 


the 4 X 4-cluster. 




k 




(0,0) 


(f,o) 


(7r,0) 


/ TT TT \ 

V2' 2) 




(7r,7r) 


n{k) 




0.9846 


0.9702 


0.9497 


0.9445 


0.9455 


0.7921 
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